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We derive anomalous equations of state for nonrelativistic 2D complex bosonic fields with contact 
interactions, using Fujikawa’s path-integral approach to anomalies and scaling arguments. In the 
process, we derive an anomalous virial theorem for such systems. The methods used are easily 
generalizable for other 2D systems, including fermionic ones, and of different spatial dimensionality, 
all of which share a classical SO(2, 1) Schrodinger symmetry. The discussion is of a more formal 
nature and is intended mainly to shed light on the structure of anomalies in 2D many-body systems. 

The anomaly corrections to the virial theorem and equation of state - pressure relationship - may be 
identified as the Tan contact term. The practicality of these ideas rests upon being able to compute 
in detail the Fujikawa Jacobian that contains the anomaly. This and other conceptual issues, as 
well as some recent developments, are discussed at the end of the paper. 
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I. INTRODUCTION 


A quantum anomaly is the loss of a symmetry in a physical system as a consequence of rigorously defining its quantum 
description. They were discovered in particle physics, where the mandatory regularization and renormalization of 
quantum field theory spoiled some existing classical symmetries M- Anomalies may also occur in nonrelativistic 
systems, but they are less well known. R. Jackiw, one of the original co-discoverers of anomalies in relativistic 
systems, very clearly pointed out the possible existence of anomalies in quantum mechanics, specifically, in systems 
with attractive inverse square potentials, or in 2D systems with contact interactions (delta function) (4j. These 
systems exhibit a classical SO (2,1) conformal symmetry, which is lost after regularization and renormalization of the 
original theory is required to formally define the quantum theory (the method of self-adjoint extensions of operators 
is an alternative way of dealing with these systems jH]. In what follows we are going to work within the regularization 
and renormalization framework). After Jackiw’s seminal paper a number of authors developed these ideas further, 
concentrating mainly on formal aspects 1 , although a nice application of these concepts was developed in molecular 
physics ||. Recently, anomalies have been investigated in ultracold, diluted trapped 2D systems. In particular, 
following the seminal work by Pitaevskii & Rosch (P&R) Q, J. Hofmann showed the possible existence of anomalies 
for 2D trapped Fermi gasses with Hamiltonian 0 
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Despite the presence of the trapping term f cP.x muJ ° x i/j] T tp a (x) in Eq. (1), which explicitly breaks the scale invariance 
of H, P&R showed that the system possesses a classical hidden SO( 2,1) symmetry, described in terms of the Lie- 
algebra generators H , D (generator of dilations) and C, with 


D = J d 2 xximji(x), 

ji = - di^ip)/2m, 

C = J d 2 x^-^ a (x). 


( 2 ) 

(3) 


Hofmann noted that due to the need for regularization and renormalization of the quantum field theory, the classical 
scaling symmetry of the system is broken, which is manifested by the existence of a ground-state energy Eb (scale) in 
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the (attractive) 2-body sector, with the ensuing alteration of the original Lie algebra. For instance, the anomaly-free 
commutation relation 


[D,H] = 2iH (4) 

is modified due to the anomaly and it now includes an “extra term” 

< 5 > 

where I is the so-called Tan contact operator EMI 

I = j d 2 xg 2 ip\ipli/jiijj t (x). (6) 

These algebraic modifications have experimental consequences that can be measured. Hofmann calculated the anoma¬ 
lous frequency shift at 8lo/ujq at zero temperature for the breathing modes and found that, in some cases, e.g., in 
the strongly interacting regime \nkpa 2 D ~ 0 ( 02 d = 2D scattering length), the change could be of the order of 10% 
compared to its symmetric classically predicted value. Other important effects due to the anomaly discussed by him 
were the shift in the ground-state energy 


E 0 = {H + uj 2 0 C) 


2u 2 (C) 


1 

47TTO 


(I), 


and in the pressure of the untrapped system (classically, due to the SO( 2,1) symmetry, £ — P = 0) 


(7) 


P = £ 


1 

4 m ’ 


( 8 ) 


where £ and I are the energy and contact density of the gas. Other similar issues and effects for ultracold Bose gasses 
have been discussed by Olshanii et al. (I3 |. and a number of papers have been addressing theoretical and experimental 
aspects of the modifications of the physics of ultracold gases due to the presence of an anomaly EM3- While much 
work has been and is being done, more work is needed in order to fully understand these anomalous effects [2lj . 
In particular, zero-temperature (anomaly) effects vs finite-temperature effects have to be further studied. Hofmann 
speculates that the experimental negative anomaly findings of the authors of [22| may be explained as a consequence of 
finite-temperature corrections. The interconnection between zero and finite-temperature effects is indeed non-trivial. 
For instance, Chafin and Schafer recently computed the grand potential H to one loop for a 2D Fermi gas at finite 
temperature (same H osc as Eq.(l), but without the trapping term) J23|. To deal with the infinities encountered they 
renormalized using the same renormalization scheme as in the description of the 2-body ground state (anomaly), and 
proceeded to calculate the shift to the free second virial coefficient to obtain 2 


Sb 2 = e pEb - 2 


dk e~ 2/3ek 



This formula exactly coincided with the well-known Beth-Uhlenbeck formula 24] 


Sb 2 = e 0Eb 




e -2/3e fe 


(9) 


( 10 ) 


2 In a theory with classical scale invariance such as this one, the phase shifts should naively all be constant, and hence the second term 
of Eqs. (9) and (10) would be absent; but the existence of scaling anomalies introduces the bound state energy scale, and things change 
dramatically. 
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< 5 (°) 

is the only non-zero phase shift for the contact potential for this problem. Eb is the finite energy of the ground 
state (the finiteness of Eb is what produces the running of the coupling constant A of Eq. (1) as a function of the 
cutoff, which in turns allows for full renormalization of Cl; that’s why one sees the presence of Eb in Eqs. (9) and 
(10)). Those equations show that there is a non-trivial comingling between zero-temperature (anomaly,Eh) and 
finite-temperature effects. It is therefore important to do a comprehensive study of regularization, renormalization 
issues in 2D systems, in both zero and finite-temperature sectors, in order to be able to interpret experimental results 
and develop a deeper understanding of these s yste ms. In this note we will propose a framework that will combine 
Fujikawa’s path-integral approach to anomalies )25ll26l| and operator methods, with an emphasis on scaling properties 
and virial theorems, in the description of finite and zero-temperature anomaly effects in 2D ultracold diluted gasses. 
In order to better understand the role of the anomaly we will consider the untrapped homogeneous system of a 2D 
Bose gas. We hope that this new framework will help us address the issues raised above. While the trapped system is 
the one that is of most interest experimentally and what motivated this work, the anomalous features are produced 
by short-distance physics and they are already present in homogenous systems as several authors have shown Ena. 
When possible, comments about trapped systems will be made, but a complete study of those systems will be left 
for future work. Just as Hofmann did using operator methods, in this paper we will also find modifications to the 
symmetric equations (virial theorem and pressure equation) due to anomalies, using Fujikawa’s approach. The extra 
term in the equations found here may, therefore, be identified with Tan’s contact term, basically Eq. (6) above. This 
will be seen in Eqs. (26), (35), and (53). However, the results in this paper are or a more formal and general nature, 
and more work needs to be done in order to fully make these identifications. The technical difficulty here is the 
explicit evaluation of functional traces of operators, which requires delicate regularization. Some progress has been 
made and it’s mentioned in the comments and conclusions. 


II. VIRIAL THEOREMS IN QUANTUM FIELD THEORY AND ANOMALOUS EQUATION OF STATE 

While quantum treatments of the classical virial theorem for particles have been available for some time, the analog 
theorems in quantum field theory are not as well known (see, however f28l | and 0)- The exception, in the case of 
nonrelativistic quantum fields (relevant for us), was the work of T. Toyoda et al. starting in the late 90s [iiOj - ljH . 
They made use of the interpretation of the “virial” operator (below) as the generator of scale transformations and 
derived a series of virial theorems for finite-temperature, nonrelativistic field theories. We will not repeat their work 
here; instead, we will present a different, path-integral derivation of their results, at least in the case of a complex 
scalar field with a quartic interaction. In their work, Toyoda et al. introduce an external potential to confine the 
system to a volume V. We will not do so, and will simply assume the system has a finite, but very large volume V 
so that the usual thermodynamics assumptions can be made. For the sake of clarity we chose this specific bosonic 
system, but the methods developed will be of a general nature for interacting fields with classical SO( 2,1) symmetry 
(the fermionic case can be treated in a similar fashion). 


A. Virial Theorem 


Consider the partition function for our system: 


Z = tr 

= J [dip*][dip]e~ s ^, 


( 11 ) 


where 

S E P = J dr J d2 x ip*(T,x) ■ VV>(r,f) 

+ g(i/j*{T,x)ip(T,x )) 2 - /j, iP*(t,x)i/>(t,x) ) . (12) 


The fields obey the usual periodicity condition in Euclidean time r. Following Toyoda et al., we’re going to demand 
that under a rescaling of the coordinates x —> Xx, the number operator be invariant: 


4 


ip* (x)ip(x)d 2 x = ip'*(x')ip'(x')d 2 x', (13) 

which determines the transformation for the field variable (since the Euclidean time will play no role in what follows 
although it will later - we will omit the time variable for the time being): 


ip(x) —> ip'(x') = A 1 ip(x). (14) 

Next, we have to make assumptions about the transformation properties of the path-integral measure under the 
scaling transformations before proceeding further: 

i) [dip*][dip] = [dip'*][dip'\ (No Anomaly) . (15) 

ii ) [dip*][dip\ = J x \dip'*][dip'] J x 1 (Anomaly) . (16) 

where J x 1 is the Jacobian for a change of variables (ip, ip*) —> (ip 1 , ip'*); J 1 = 1; see appendix. 

Consider case i) first. Make the transformation ip(x) = Xip'(Xx). By relabeling the integration variables, it’s easy to 
see that only two terms change in the exponential in the path integral: 


1 

2 


d 2 x Vip*(x) ■Vip(x) 


2 


d 2 x'Vip'*(x')-Vip'(x'), 


(17) 


g J d 2 x (ip* (x)ip(x)f = gX 2 J d 2 x' (ip'* (x')ip'(x')f . (18) 

Make an infinitesimal scale transformation A = 1 + 77 . Then the partition function becomes (the r dependence has 
now been restored due to the r integration): 


J [diP*][diP]e- SE ^^ = j[diP'*}[diP']e- < < SE ^^ +ss ) 

= j [dip'*][dip']e~ SB ^ > x jl ~ j J d 2 x , 'V'ip'*( T ,x') ■'V'ip'( T ,x') 
- 2gg J dr J d 2 x' (ip’*(r, x')ip'(T, x ')) 2 | . 

Therefore, in the large volume limit 3 


(19) 



d 2 x -Vip^(r, x) ■ Vip(r, x) 




d 2 x g(ip\r,x)ip(T,x)) 



( 20 ) 


where 


(A) = Z- 1 tr (e-^n-^A) . 


(21) 


Now, using 


3 At this point there is no need to make a distinction between primed and unprimed variables. 
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i/j(t,x) = e T ^i/j(x)e , 

ip\r,x) = e T ^il>^ (x)e ~ T ^, (22) 

K = H- iiN, 

we get 

(Vif>^(r,x) ■ S7 'i/j(t,x)') = Z~ x tr e T ^X/^ (x ) ■ Vilj(x)e ~ T ^ 

= Z- 1 tr (e~ Tk e- 0k e Tk y^(x) ■ Vtp(x )) ( 23 ) 

= (V?/v (x) ■ Vtp(x)) : r-independent. 

A similar treatment shows that the RHS of Eq. (l20l) is r-independent, which produces an overall /3 factor on both 
sides, hence finally giving 


/ ^Vip^x) ' W’Ozy = - (\j d 2 xg{^(x)tlj{x)) 2 '^ , (24) 

or ( H 0 ) = (K) = —(H int ) = —(V): Virial Theorem. (25) 

(Recall that in the first-quantized version, 2 (K) = n(V), for V = ar n or for V(Xr) = X n V(r), such as 8 2 {r) ( n = —2).) 

In case ii), we have to add the contribution due to the Jacobian, Eq. m- Following similar steps to the previous 
section, one obtains a modified virial theorem (see appendix): 


(K) 


A 

— ( V ) + — : Anomalous Virial Theorem, 


(26a) 


where (Tr is a functional trace) A = —Tr (1 + x ■ V)J(r — t)S 2 (x — x) 


is the anomalous contribution from the 


Jacobian for the infinitesimal transformation ip(x) —»• ip'(x) = ip{x) + 5^(x), S^(x) = —r]( 1 + x ■ \7)ip'(x). 


Although the emphasis in this paper is the study of homogeneous systems, the derivation of the anomalous virial 
theorem including a trapping term of the form / d 2 x x 2 ^( r , x x ) along these lines would readily give 
the following version: 4 


(K) = -(E) + J + (J ^Y-dfx^tx)^)^ . (26b) 

Comparison with Eq. (7) would suggest that the anomaly term A should be identified with Tan’s contact in the 2 D 
case (see also Eq. (2) in ref. 0 )- 

Evidently, an ex plic it evaluation of A would require careful regularization (and possible renormalization) of the 
functional trace ([33(; see below). 


B. Anomalous Equation of State 

Consider the partition function for the scaled system (x —> x' = Xx , r —> r, i/j(t,x) —>• ip'(T,x r ) = A -1 , 0(t, x)) 


4 Notice the thermal expectation values in Eq. (26b) are computed using the full Hamiltonian of Eq. (1). 
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Z x = J [dtp *} , [dilj} , e- s ^ f, ' X 

where 


(27) 


Se P ' X = J dr J d 2 x’ U • V>'*(T,f') 

+g (lp'(T, x')lp'*{T, X 1 )) 2 - /X 1p'(T, x')lp'*{T, X ') ) 

= J dr J d 2 x ^tp*{T,x)-^-tp{T,x) + A~ 2 (^\/tp*{r,x) ■ V^(r,x) 

+5 W>* (u x)tp{r, x )) 2 ) - ntp* (r, x)tp(r, x) J. 


(28) 


Using ((11) and (16)) we get 

Z A = J A tr ^ e -^( A “ 2ff - /iiv )) . (29) 

(H = H 0 + H int = i J d 2 xWtp^{x) ■ Vtp(x) +g J d 2 x (ip^(x)ip(x)) 2 ) 

Z x must correspond to a change of volume in the system (not in temperature or pressure) 5 . In the large volume limit 


Z x = e~^ X = e d pA \ 

A = 2D “Volume” ; A a = scaled 2D “Volume” = A 2 A. 
For infinitesimal scalings A = 1 + 77 , we get (AA = 2r/A) from Eq. (29) 

Z 1+ v = e P p ( A + AA ) _ j 1+, ?tr . 

As before J l+V gives an infinitesimal contribution (see appendix) 

J 1+r > = 1 + 2 r/Tr (§ s 6 3 {x - y)) , 


(30) 


(31) 


(32) 


and 


0 s = -(l + x-V). (33) 

For a 2D volume A, under x —> Ax, A —>• A +AA, AA = 2?yA. Expanding Eq. (31) in 2g on both sides, and collecting 
terms we get Toyoda’s result with the anomalous contribution: 


PA — {Ho) + (Dint) + 


Tr ( 6 s 5 3 {x-y) S j 

J3 


(H) + 


Tr ( 6 s 5 3 (x-y )) 

0 


(34) 


5 See [3^1 for a graphical proof of this statement. 





7 


which can be rewritten as 6 


Tr 


E-P= — 


(e s 5 3 (, 


x-y 


PA 


£ = 


(£) 

A ' 


(35) 


Equation (35) clearly displays the anomaly effect on the thermal (many-body) system: if the symmetry is preserved 
(J = !),£ = P, as is expected. But in general, if Tr ( 9 S 5 3 (x — y )) ^ 0, the equation is modified 7 . 


Again, comparison with Eq. (8) and Eq. (3) of ref. HU suggests as before that the term 
identified as Tan’s constact term in the 2 D case. 8 


Tr {S s S a {x-y)) 
13 


is to be 


C. Algebraic considerations 

Toyoda et al. used a mixture of scaling arguments and algebraic manipulations to deduce their equations of state. 
We will now explore the algebraic content of our formulation, and compare with theirs, as well as with Hofmann’s. 
Using their confining potential method, they arrived at the following result (their formula (2.5) of reference [61] and 
the similar one for the 3D case in a previous paper differ from ours by an overall factor. They use isotropy arguments 
to derive their equations of state, which agrees with ours. We don’t need to make any extra assumptions): 

PA = ~([G,H]) = (H), H = Ho + H int . (36) 

G is defined in Eq. (39) below. We can obtain the same formula within our formulation (in the non-anomalous case): 


Z 1+T l = tr e -PPH+SH)-v,N) = tr e -P(H-iJ.N) ^ _ p SH j ( 6H ^ H jV] = 0) 

= Z 1 -1 3 tr (37) 

= Z 1 + Z\2r]p)AP. 


Therefore, 


But, 


AP=~{SH). 

2 y 


(38) 


SH = H 1+ri - H 1 = -2 yH = [G, H]rj , 

G = “virial” operator (generator of scalings) = — 


d 2 x x ■ ijp ( V — V ) tf>, 


(39) 


6 In thermal equilibrium and in the large V (A) limit, one expects Tr(6J s 5 3 (a; — y)) to be independent of r and x in our case (untrapped), 
which will provide an additional j3A factor that will cancel the similar term in Eq. (35). This indeed happens explicitly in the Jacobian 

calculation of paper |33| . If one includes the trapping term f d?x —^one works with the anomalous density of Eq. (35). 

7 While we have used the terminology of quantum anomalies, strictly speaking, since Eq. (13) is not necessarily a symmetry of the action, 
the formal Eq. (35) should be considered more of a “modified” equation vis-a-vis Toyota el al. Furthermore, since the treatment in both 
this and their work is basically formal, after proper regularization and renormalization (presumably somewhat different in each case), 
both might still give the same physical results, despite the difference in the formal equations [3^1 . In the case of an actual SO( 2,1) 
quantum anomaly (section D), the presence of a non-trivial Jacobian J / 1 would signify a real physical effect (symmetry breaking). 
The path-integral methodology of the current approach is essentially identical in both cases though, and hence the early use of the 
language of anomalies. 

8 This method also gives the results for D = 3. 






and we have used [G, H] = —2 H, which can be explicitly verified for our type of Hamiltonian, and using the standard 
commutation relations for ip(x), ip^(x). 

Using Eq. (1551) and Eq. ( 1551 ) we then arrive at Eq. m- This is a self-consistent derivation, if one ignores the 
zero-temperature 50(2,1) symmetry of the system and its possible anomaly, but as we saw with the example of the 
calculation of the second virial coefficient of the 2D fermion gas, only at one’s peril! Indeed, if we do assume that 
our virial G is part of the set of 3 generators of 50(2,1) that generates dilations, one has to take into account the 
anomaly in deriving consequences of the algebra (as Hofmann also showed). The full dilation operator D is related 
to G as follows: 


D = -2tH - iG. (40) 

In the non-anomalous case, the commutator with H is known to be (as part of the 50(2,1) symmetry, where both 
time and space are transformed) 


[D,H]=2iH. (41) 

Eq. (I4T1) is obviously equivalent to [G, H] = —2H. So Eq. ( 1561 ) can be written as 


PA = - l -([D,H]). (42) 

In the anomalous case, [ D , H] = 2iH + extra terms , which in Hofmann’s case of the 2D fermion system gave the 
corrected equation, Eq. (8), and we of course also anticipate this in the bosonic case described here. Embedding the 
scaling operator G into the full 50(2,1) via Eq. (l40l) is an expedited way to find a formal corrected Eq. for the rela¬ 
tionship between pressure and energy density, but it’s not too insightful in computational terms, which we would like 
to have in order to take such equation into an actual equation of state, virial expansion, etc. Therefore, a path-integral 
treatment similar to the one described above, but now performing both time and space scalings, t —> A 2 t, x —> Ax, 
should be the way to set up the path-integral derivation of the anomaly-corrected equations. We now turn to this task. 


D. Derivation of 50(2,1) anomalous equation of state 

In the previous subsection we considered scalings only in the spatial coordinates in order to make contact with Toyoda 
et al.’s derivation of their equation of state (what they called virial theorem), and we already saw that on account 
of the scale anomaly their equation had to be modified by the addition of a Jacobian piece. In order to display the 
anomaly of the full 50(2,1) symmetry we need to consider the partition function with both time and space scaled 
accordingly, with the ensuing field transformation 9 : 


Z A = J[dil)'*}[dil>']e- s E ' 5 Vb’/’ ,] ) ( 43 ) 

where the primed (scaled) variables are: 

r —>• t' = A 2 t, (44) 

x —» x' = Ax, 

ip(r, x) — > ip'( T ', x 1 ) = &)■ 


Infinitesimally, for A = 1 + rj 


9 


Notice the field transformation remains the same as in the case when only the spatial components are scaled; this may have significant 
relevance in regularizing the Jacobian [33ll . 
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St = 2 ?yr, 

Sx = r/x, 

Sip = —r] (l + x ■ V + 2tS t ^ ip(r, x). 

Performing explicitly the scalings in the action gives (see Eq. (12) and notice that P' = A 2 /3) 

si' f) 'W*A , ] = s x z^WM- 

Z x can therefore be expressed as 


(45) 


(46) 


Z A = J x J [dip*][dip] e 
= J A tr ( 'e -^ H ~ A2 ^)) . 


-si 


In the large volume limit, Z x becomes 


(47) 


Z x = e f) ' p ' A ', (48) 

where now, due to the additional rescaling of t, we expect changes, not only in the volume A, but also in /3 and P. 
These changes are contained in both the Jacobian and the trace factors of the RHS of Eq. (47). 

We now make an infinitesimal change A = 1 + rj and compute the changes on the left and right hand sides of Eq(47). 
On the LHS: 


SZ 


e p ' p'a' - e 0PA = e^ FA (pPSA + PAS/3 + P^-A8P) 

op 

dP 

= Z(4r)f3PA + 2r)/3—A/3). 


(49) 


From PA = ^lnZ we get PA^ = —PA + and hence, 


SZ = Z(2 V P) ^PA + = Z(2r 1 p)(PA - ( H) + »(N)), 


where we used = —(H) + /J,(N). 

From the right hand side of Eq. (47) we get 


(50) 


SZ = Z (2r]Pn(N) + 2rjTr (§S 3 (x - y) 
where now, in view of the transformation Sip, Eq.(45), the anomaly operator includes a time piece 

9 = - (l + 2 rd T + x ■ v) . 

Therefore, the SO( 2,1) anomalous equation is (as before £ = ^-) 


(51) 


(52) 


Tr 


£-P = — 


(^ 3 ( 


x-y 


PA 

Eqs. (26), (35) and (53) are the main results of this paper. Comparison with Eq. (8) and Eq. (3) of ref. 
suggest that in 2D, should be identified with Tan’s contact term. 


(53) 

would 




10 


III. COMMENTS AND CONCLUSIONS 

We have derived in this paper anomalous-modified equations of state for 2D scale invariant systems with contact 
interactions using scaling arguments and Fujikawa’s path-integral approach to anomalies, Eqs. (35) and (53). In the 
process, as a by-product, we also obtained an anomaly-corrected virial theorem for such systems, Eq. ( 26). The work 
here is of a formal nature, and it is meant to illustrate how in principle one can account for the anomalous effects in 
a clear fashion, which we hope will help us address some of the issues discussed in the introduction. Indeed, work 
based on this paper to investigate those topics is in progress, and it has already produced some concrete results (M-, 
see below). A few comments are in order: 

• While we chose the relevant 2D contact term interaction for bosonic fields as a concrete case in this paper, most 
of the calculations and formal manipulations rely only on the scale transformation properties of the two-body 
potential, V(Xr) = A ~ 2 V{r). Therefore, the equations derived here also apply in the case of inverse square 
potentials in arbitrary spatial dimension and other SO{ 2,1) invariant systems such as anyons [37l - l42| . 

• Likewise, the final results for the Eq. of state and the virial theorem apply to fermionic systems. 

• The emphasis in this paper was on the impact of anomalies in 2D, or in general, systems with a classical SO( 2,1) 
symmetry; however, most of the calculations and the ideas developed in our framework apply even when there 
is no such classical symmetry, or when the symmetry exists at the quantum level in some asymptotic regime (as 
in atoms at unitarity 0). In particular, our approach produces Eqs. similar to Eqs. (35) and (53) in D = 3, 
agreeing with Toyoda et al.’s results, if one ignores the Jacobian term (see footnote 7). This procedure also 
seems to produce Tan’s contaact term in D = 3. 

• As it has been stated above, comparison of Eqs. (26), (35) and (53) with the literature would suggest the 
identification of the Jacobian contribution with Tan’s contact term, and in the 2D case, its presence signaling 
the existence of an 50(2,1) anomaly. In 2D, it appears as though one need use only the spatial scaling 
transformation to derive the contact term, regardless of the existence of a classical 50(2,1) symmetry. When 
dealing with the issue of a potential 50(2,1) anomaly, one has to consider both time and space dilations. The 
operators 6 S and 8 so obtained are different, but a proper treatment of the functional traces should produce the 
same answer. Further study needs to be given to this issue. 

• It is clear that the practical relevance of the approach presented here rests upon one’s ability to calculate the 
Jacobian, i.e, to carefully regularize it to produce usable finite (possible bare) answers. We have made progress 
in this direction in the bosonic sector, for the homogeneous case discussed explicitly in this paper (without the 
harmonic trapping term) jUj, and we have made contact with previous results in the literature that confirms 
the validity of our approach in this case 0, [45[ . We now give a brief summary of the salient features of the 
results in 133! : Following the path-integral version of Noether’s theorem, and using the scaling transformations, 
Eq. (A.2) of the appendix, one arrives at the expression containing the anomaly term on the right hand side: 


<V) = -*Tr (85\x-y)I 2 ), (54) 

with 6 = 1 — x • V — 2 tdt'j , I 2 = ^, and Tr refers to both matrix and functional traces. 

In the case of a constant background (consistent with the homogenous case), we can ignore the derivative 
contributions in 8 (notice we’re using the notation x = ( xo,x )) and regularize Eq. (54) using the eigenbasis 4> n 
of a Hermitian operator M as follows: 


S 3 (x - y)I 2 = 'Y^(t> n {x 0 ,x)(t)i{y 0 ,y) -A 5\{x - y)I 2 = (-p) (f> n (xo,x)4>\{yo,y), (55) 

n n ' ' 

with the property that i?(0) = 1 so that at the end of the calculation we send A —> 00 and lim R (-M-) = 1. 

A—>oo '’ 

The trace is then regularized as 
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Tr [5 3 r (x - y)h] = Tr R S 3 (x - y)I 2 

with the choice of the Hermitian M determined by the one-loop structure of the theory 


(56) 


M _ lid t + \+g- 2gip*ip + ie ^ -g^ 2 \ 

\ -g4>* 2 —id t + ^ + /j, - 2g'i/j*i/j + ie) ’ 


where the fields in M are constant background fields. A class of regulators of the form 


(57) 


R 




(58) 


gives, for the zero temperature case 


(d,r) 


g 2 (Vi O 2 

4-7T 


and, in the finite-temperature case considered in this paper, Eq. (53): 


(59) 


2£ — 2P = — 


9 2 {V^) 2 

47T 


(60) 


In both cases we find agreement with the literature [3, More work remains to be done, including the case 
when the number density is not constant, in order to make contact with other realistic situations of interest in 
atomic and molecular physics, as well as in condensed matter physics, and possibly in applications to holography 
and the AdS/CFT duality 
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Appendix: Jacobian contribution to Eqs. (26), (35) and (53). 

In the change of variables, Eq. (16), we need to calculate the Jacobian 


( 5ip(x) 8'ip(x) \ 

115 ggg Udet (S(x,y)), (A.l) 

<5b'(y) <5 b'*(y)/ 

for infinitesimal scalings A = 1 + g. We are using the notation x = (r,x),y = (f, y); S(x , y) is the 2x2 matrix defined 
by Eq. (A.l). 


The general 5*0(2,1) infinitesimal dilation is 
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5x = rjx, 

St = 2 ryr, 

Sip = y9ip(r, x), 
Sip* = y9ip*{r, x), 


(A.2) 


where 


0 = -(l + £- V + 2rd T . 


(A.3) 

For subsection II-A, we need consider only the spatial part, 9 S = —(1 + x ■ V). For these transformations S(x,y) 
becomes 


5 3 {x - y) - r]9 s S 3 (x - y) 0 

0 S 3 (x-y) - r]9 s 5 3 (x - y) 

5 3 {x-y) 0 \ f-r]9 s S 3 (x-y) 0 

0 5 3 {x-y)) + \ 0 - r)9 s 5 3 {x-y ) 

= l + r]B(x,y). 


In Eq.(A.4) S 3 {x-y) = S(t - r’)S 2 (x - y), B = ^ ® sS ^ ^ s^x - y)) ' 

We now use the well-known identity (Tr includes a functional trace as well as the 2x2 matrix trace) 


where Tr is the functional trace only. 
Therefore, for A = 1 + y (A -1 = 1 — 77 ) 


det S = e 


Tr In S 


= e Trln(I +riB) 

= 1 + 2t]Tt{B) 

= 1 - 2r?Tr (§ s S 3 (x - y)^j 


(A.4) 


(A.5) 


J 1 -" = 1 - 2r]Tr ( 9 s 5 3 {x - y)) . 


(A. 6 ) 


For the change of variables in the opposite direction, primed to unprimed (used in subsections II-B and II-D), we 
need J A , so the trace term in Eq. (A. 6 ) is positive in those cases. In subsection II-D, we need to replace 9 S by 9. 
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